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Abstract. - Network theory provides a rich toolbox consisting of methods, measures, and models 
for studying the structure and dynamics of complex systems found in nature, society, or technol- 
ogy. Recently, it has been pointed out that many real-world complex systems are more adequately 
mapped by networks of interacting or interdependent networks, e.g., a power grid showing inter- 
dependency with a communication network. Additionally, in many real-world situations it is 
reasonable to include node weights into complex network statistics to reflect the varying size or 
importance of subsystems that are represented by nodes in the network of interest. E.g., nodes can 
represent vastly different surface area in climate networks, volume in brain networks or economic 
capacity in trade networks. In this letter, combining both ideas, we derive a novel class of statis- 
tical measures for analysing the structure of networks of interacting networks with heterogeneous 
node weights. Using a prototypical spatial network model, we show that the newly introduced 
node- weighted interacting network measures indeed provide an improved representation of the un- 
derlying system's properties as compared to their unweighted analogues. We apply our method to 
study the complex network structure of cross-boundary trade between European Union (EU) and 
non-EU countries finding that it provides important information on trade balance and economic 
robustness. 



Introduction. — Complex network theory has been 
. shown to be a powerful tool for analysing the structure and 
• • function of many complex systems in nature, society, and 
. ^ technology. Various kinds of measures have been defined 
recently, mostly based on counting nodes, paths, links or 
5— (triangles in a network [l}|3]. The field of application is 
^wide-spread, e.g., considering social [3], trade [H], biologi- 
cal (6[|7j, communication [s] and climate networks [9 13 



When applying network theory to real- world networks, it 
is often not sufficient to describe the underlying complex 
system by an isolated network. Instead, a network of inter- 
acting networks may provide an improved representation 
as was shown, e.g., in an analysis of the interdependency 
between the Internet network and the Italian power grid 
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during a blackout in 2008 14 . In general, interdependent 
networks behave much differently from isolated ones in 
terms of robustness to random failure, expected network 
properties 15 16 as well as synchronisation behaviour 
17 . In order to quantify the structure of interdependent 
networks, the recently introduced interacting networks ap- 
proach compromises a set of cross-network measures in 
analogy to the canonical network measures designed for 
isolated networks 12 18 19 . The method has been ap- 



plied successfully for analysing the dynamical structure of 
the lower atmosphere by constructing coupled climate net- 
works from pairs of geopotential height fields at different 
isobaric surfaces 



12 18 



When applying network theory to real- world problems, 
nodes need not all bear the same importance for the net- 
work's properties, but there may be nodes having a strong 
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Fig. 1: (Colour online) A set of three subnetworks Gi — 
{Vi,Ei,) (green), G, = {V,,Ejj) (blue) and Gk = iVk,Ekk) 
(orange) connected via both internal (solid lines) and cross- 
(dashed lines) links. Each node v represents a share of the 
whole domain of interest and is therefore provided with an in- 
dividual weight Wv The node sets Vi and Vj share a common 
boundary in some generalised space, while Vk remains sepa- 
rated. 



impact on the network's topology that is not reflected by 
link properties alone. E.g., nodes can represent vastly 
different surface area in climate networks [9 



11 



13 , vol- 



ume in brain networks [6j[7| or economic capacity in trade 
networks I^O^jST]. Therefore standard network measures 
treating all nodes the same may not represent appropri- 
ately all structural properties of the underlying complex 
system (in the following referred to as the domain of inter- 
est) [22j. To take this into account, it is useful to assign an 
individual weight to every node according to the share of 
the whole domain of interest that is represented by that 
node. The concept of node splitting invariance (n.s.i.) 
yields a set of measures allowing for a precise estimation 
of network characteristics when dealing with inhomoge- 



neous node weights 13 



The main aim of this work is to combine the frameworks 
of n.s.i. and interacting network measures in order to pro- 
vide a general tool for investigating interaction structure 
within a network of networks with heterogeneous node 
weights (Fig [T]). Starting from the cross-network mea- 
sures |12| , we introduce a general construction mechanism 
to derive n.s.i. cross-network measures. Subsequently we 
apply these measures to a global model network with a 
structure prototypical of spatially embedded networks of 
networks and validate the results for n.s.i. cross-degree k^* 
and n.s.i. cross-clustering coefficient CI* . To demonstrate 
the potential of our approach, we carry out an analysis 
of a trade network divided into two subnetworks, EU and 
non-EU countries, where nodes represent individual coun- 
tries weighted according to their Gross Domestic Product 
(GDP) as a proxy of economic power. 

Methods. — We consider a network G — {V, E) with 
a given set of nodes V , links E and the number of nodes 
= As we identify every node v € V with a nat- 
ural number p < N, the network G is represented by its 



adjacency matrix A with Apq = 1 if {p,q} G E, Apq = 
Oif {p,q}^E. 

Let this network be divided into M > 2 subnetworks 
Gi = (yi,Eii),i = 1,2,...,M with subsets of nodes Vi, 
such that V = IJ,; with n V,- = for i ^ j. The set 
of links E then splits into two types of links. There are 
sets of internal links En that connect nodes v d Vi within 
a subnetwork Gi and cross-links Eij that connect nodes 
V E Vi with nodes q G Vj in the subnetworks Gi and Gj , so 
that E = [J- J Eij (Fig. [l|. Notice that a network might 
be divided into several subnetworks by splitting a con- 
tinuous domain of interest into smaller sub-domains, e.g., 
when dividing a landmass into adjacent countries forming 
subnetworks {Gi and Gj in Fig.[T]). However, it might also 
occur that the two sub-domains are separated, since they 
display different observables or layers, such as a power 
grid with power plants as nodes and transmission lines as 
links which show interdependency with a communication 
network, where servers are represented by nodes and links 
indicate wires between computers {Gk is separated from 
G, and Gj in Fig.[l| [m]. 

Assuming that every node v € V represents a part of 
a larger (maybe continuous [23| ) domain of interest, we 
assign it a weight Wy > which represents its share of 
the domain of interest. Therefore, when redefining any 
unweighted network measure, nodes should no longer be 
regarded as having a well-defined location, but their pos- 
sible locations are constrained by the boundary of the rep- 
resented area Wy. This implies that hypothetically split- 
ting one node v £ Vi into two adjacent and similar nodes 
v' € Vi and v" £ Vi the former node weight Wy should split 
up into Wy' and Wy" such that Wy — Wyi -\- Wy" . Since we 
assume links in a network to represent similarity or asso- 
ciation between pairs of nodes, we treat v' and v" as being 
connected to all neighbours of the former node v and to 
each other. Given that the sampling G already provides 
a good approximation of the underlying complex system's 
structure, this node-splitting operation merely refines the 
network representation and should not affect the network 
measures much. Taking this into account yields a four- 
step construction mechanism for transforming unweighted 
network measures into their weighted counterparts as sug- 
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(a) Sum up weights Wy whenever the unweighted measure 
counts nodes. 

(b) Treat every node v G V a.s connected with itself. 

(c) Allow equality for v and q wherever the original mea- 
sure involves a sum over distinct nodes v and q. 

(d) "Plug in" n.s.i. versions of measures wherever they 
are used in the definition of other measures. 

(b) and (c) are derived from the fact that when splitting 
one node into two similar ones, the new nodes will be 
connected because of their similarity, (d) can be seen as 
an inductive step when applying (a) - (c) to an unweighted 
network measure. Treating every node as connected with 
itself, we additionally introduce the extended adjacency 



matrix A^^ — Apq 



Spq, where Spq is Kronecker's delta. 
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n.s.i. cross-network measures. — Here, we derive 
two informative measures for analysing networks of net- 
works with heterogeneous node weights: (i) n.s.i. cross- 
degree and (ii) local n.s.i. cross-clustering coefficient. Fur- 
ther metrics can be easily derived analogously (see [13| ). 

n.s.i. cross-degree. The unweighted cross-degree [12| 



(1) 



measures the number of nodes q ^ Vj that are connected 
to a node v €z Vi (typically i ^ j). Here we use mechanisms 
(a) and (b) to construct the n.s.i. cross-degree as 



ki* = J2^' 



(2) 



The n.s.i. cross-degree no longer just counts the number 
of nodes that v is connected to. It rather measures the 
share of the whole sub-domain of interest given by the 
subnetwork Gj that the node v G Vi'is connected to. While 
the standard cross-degree can only take integer values in 
the range of 0,...,N — 1, the n.s.i. cross-degree fcj* can 
assume real numbers in the range of < kl* < Wj , where 
Wj = J2qev total weight of all nodes qinVj, e.g., 

the area on the whole sub-domain of interest represented 
by the network Gj (Fig. [T| . 

n.s.i. local cross- clustering coefficient. The un- 



weighted local cross-clustering coefficient 12 



(3) 



quantifies the probability that two randomly drawn neigh- 
bours p,q £ Vj oi V e Vi are also linked. Illustrating 
the construction mechanism, we convert the local cross- 
clustering coefficient into a n.s.i. cross-measure step by 
step: 



ki{ki 1) p^q^y. 

^ ^tpWpA+qWqA^ 
p,qGVj 



[kiY 



(d) 1 



{ki*Y 



J2 AtpWpAj^WqA+, = cr e [0,1]. (4) 

p,qeVj 



This only works if kl* > 0. Otherwise one may consider 
Gi* as undefined or set Gi* = 0. The n.s.i. local cross- 
clustering coefficient gives the probability that two ran- 
domly chosen points on the continuous domain of interest 
of the subnetwork Gj which are neighbours of a point in 
the area represented by the node v & Vi are also neigh- 
bours. 




o o o o 




Fig. 2; (Colour online) Model of a spatial network G = {V, E) 
consisting of two subnetworks Gi = (yi,-Eii) (filled circles) 
and G2 = (V2,-E22) (open circles) representing layers parallel 
to the Earth's surface and links between nodes put randomly 
according to the distance between them. Only cross-links (red 
solid lines) emerging from two nodes in the upper layer are 
shown. 



Application 1: Spatial network model. — As- 
sume a global model network G = {V, E) consisting of 
two subnetworks Gi = {Vi.En) and G2 = (V2,-E22) rep- 
resenting layers parallel to the Earth's surface with a dis- 
tance d ^ R, where R is the Earth's mean radius, as a 
typical representative of coupled climate networks 
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and other spatial networks [24j (Fig. [2]) . The node distri 
bution in each layer is chosen to be identical and nodes 
are distributed globally with angular distance of 2.5° in 
latitude and 5.0° in longitude yielding a total number of 
A'' = 10, 658 nodes. Links between nodes are considered 
to display a certain degree of similarity meaning that ge- 
ographically close nodes are more likely to be linked than 
those spanning a larger distance. However, a small num- 
ber of long-distance links such as teleconnections in cli- 
mate networks [9] or intercontinental flights in airline net- 
works [25| may also be present. Therefore, in our model, 
links between all pairs of nodes v and q are introduced 
with a probability P^s^q) o: exp(— s^^/A), where s-uq is 
the geodesic distance between v and q. The typical length 
scale A = 1000 km was chosen to yield a link density of 
p « 0.05. Regions close to the poles are expected to dis- 
play larger local link density due to the increased node 
density in these areas (Fig. [2]). The given setting pro- 
vides an idealisation of typical spatially embedded net- 
works 24 , where the nodes are located on a geographical 



grid having not necessarily a homogeneous node density. 
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Fig. 3: (Colour online) Zonal averages of local standard cross- 
network measures and the corresponding n.s.i. cross-network 
measures for a spatial network model. (A) N.s.i. cross-degree 
k^*{9) and (B) n.s.i. local cross-clustering coefficient C^*{9) 
provide an improved representation of the underling model 
system's properties when approaching the poles {9 — >■ ±90°). 
Symbols represent means, error bars standard deviations based 
on an ensemble of f 00 realisations of the network model. 



In this case, nodes v at different latitudes 9^ represent a 
different share of the Earth's surface. We assign an indi- 
vidual weight Wy — cos{9y) to every node v so that nodes 
close to the poles have less weight than nodes near the 
equator. This choice of weight approximates the area ac- 
tually represented by a node in an angularly regular spher- 
ical grid [9 13 . 



In analogy to studies showing that the geometrical 
structure of a spatially embedded network strongly in- 
fluences the distribution of network measures 26 27 , we 
expect that local network properties, which are mainly 
based on counting adjacent nodes, should show remark- 
ably increased values in regions with high node density. 
We compute zonal averages for the standard cross-degree 
k'^{0) and the n.s.i. cross-degree k1*{9) for all nodes at 
latitude 9 (only the results for Gi are shown due to sym- 
metry). The values of kfj{9) rapidly increase towards the 



poles (Fig.[3j^), while due to the inhomogeneous node den- 
sity, values are low at the equator. In comparison, k'^*{9) 
which measures the area on the surface of G2 that is con- 
nected to V exhibits a constant value even close to the poles 
(Fig. [Sj^). Here, the increased node density in this area 
compensates with the comparatively small node weight. 
Since we have considered a homogeneous network with 
no nodes having any outstanding properties, but pairs of 
nodes being randomly connected according to their dis- 
tance, kl* (9) represents the network's topology much bet- 
ter than kl{9). 

The behaviour of the standard local cross-clustering co- 
efflcient Cy{9) is similar to that of kl{9) and the values of 
Cy{9) increase rapidly towards the poles (Fig. |3|3). Due 
to the increasing node and local link density at the poles, 
nodes in this area have far more neighbours in their closer 
geographical surrounding than nodes near the equator. 
Since those neighbours are also more likely to be found 
close to the pole, it is quite probable that two neighbours 
of a node near a pole are also neighbours, yielding higher 
values of Cy{9). Cy*{9) is more robust to these circum- 
stances showing an almost constant distribution of values 
over the whole globe. However, Cy{9) and C^*{9) do not 
converge in the equatorial region. This disparity can be 
dealt with by further adjusting the n.s.i. cross-network 
measures and introducing a typical weight to in the defi- 
nition of the measures 
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Application 2: International trade network. 

Trade relationships between actors in world trade have 
recently been intensively studied 20pT 28 and, in a glob- 
alised world, provide a particularly interesting example of 
a complex spatial network of interacting networks. The 
following analysis of a trade network illustrates potential 
for applications of interacting network measures pointing 
out the importance of taking into account properly cho- 
sen node weights Wy . In contrast to the above application, 
where node weights reflected the surface area represented 
by nodes, here we are interested in node weights indicat- 
ing the national economic power that is involved in trade 
relationships. 

In particular, we perform an analysis of trading struc- 
ture across the border of the common market of the Eu- 
ropean Union (EU) by studying a network of networks 
consisting of two subnetworks: Gi representing all 27 
EU countries and G2 containing 115 non-EU countries 
(Fig. [4]) . Nodes v represent countries and links indicate a 
significant amount of trade between countries (as in [5]). 
Specifically, a link is created if the total value of the mu- 
tual reported trade in 2009 according to comtrade.un.org 
between two countries accounts for at least 5 % of the total 
trade of any of the two countries. Visualising this trade 
network indicates that defining EU and non-EU countries 
as two subnetworks is justified from the data by the fact 
that the nodes representing EU countries tend to form a 
tight cluster [l2] (Fig.|4]). Additionally, we assign a weight 
Wy to every node according to the 2008 Gross Domestic 
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Fig. 4: (Colour online) Visualisation of a trade network con- 
sisting of 142 countries that divides into two subnetworks, one 
representing the 27 members of the EU (blue) and one rep- 
resenting all non-EU countries (orange) using the force atlas 
2 algorithm in Gephi 29 . Node size is proportional to node 



weight Wv- the countries' GDP in 2008. 

Product (GDP in USD). The GDP data is provided by the 
International Monetary Fund (IMF, http:/ /www.imf.org). 
By doing so we do not only consider the total number of 
countries participating in trading but also the strength of 
different national economies. 

We compute the n.s.i. cross-degree k^* as well as the 
n.s.i. cross-clustering coefficient C^* for v G Vi and vice 
versa (Fig.jsj). The dashed line indicates the expected val- 
ues if Wv would not correct the examined network mea- 
sure at all. For the cross-degree as well as the local 
cross-clustering coefficient, we observe no linear depen- 
dency between the weighted and the unweighted network 
measures. However, for almost every node the values of 
n.s.i. cross-network measures are significantly larger than 
their unweighted analogues. This observation indicates 
that the unweighted measures underestimate the corre- 
sponding properties of the underlying trade structure (the 
domain of interest). 

The cross-degree kl counts the number of significant 
trading partners of a country v across the border of the 
common market. In contrast, fc^* takes into account each 
country's GDP, meaning its value displays the total eco- 
nomic strength that v is trading with (Fig. [5]^.). This 
way we identify the two countries with the highest devia- 
tion from the equality kl* /Wj = kl/Nj as Great Britain 
(GBR) for the EU countries and Brazil (BRA) for the 
non-EU countries. For both countries, k^* /Wj is much 
larger than k^ /Nj , implying that both have comparatively 
few trading partners across the border of the EU common 
market, but the GDP of those partners sums up to a large 
economic strength that BRA and GBR are trading with. 
This can indicate both, an either positive or negative trade 



0.6 



0.4 



0.2 







[bra] 






'0^ . .. 

<>• CQO' 


/ 

/ 

/ 

/ 

/ 


2§ / 




Y / 

1 A 


non EU countries 
oO- EU countries 



0.2 



0.4 0.6 

kl/Nj 



1.0 



B 



1.00 



0.95 



cj- 0.90 



0.85 



,0 



^4 • non EU countries 
• •• EU countries 



0.2 



0.4 0.6 

ci 



1.0 



Fig. 5: (Colour online) Scatter plots of (A) normalised n.s.i. 
cross-degree kl* /Wj and (B) n.s.i. local cross-clustering coeffi- 
cient Ci* with their unweighted counterparts for the network 
of significant trade relations in 2009. Ci* is shown for all nodes 
with kl > 1. Symbol size is proportional to node weight w„: 
the countries' GDP in 2008. For reference, dashed lines indi- 
cate equality of normalised weighted and unweighted measures. 



balance, of the country with its partners across the EU 
boundary. BRA exports comparatively many goods to the 
EU countries with strong economies, whereas GBR mainly 
imports goods from strong economies across the boundary 
of the common market. In order to obtain a broad insight 
into the trading behaviour of a country it is therefore nec- 
essary to evaluate both versions of cross-degree as they 
may be interpreted as different properties. 

The local cross-clustering coefficient gives the prob- 
ability for a country v G Vi to form a triangle with coun- 
tries p,q € Vj. In terms of trading this implies that a 
considered economy v £ Vi has, besides trading directly 
with an economy g G 1/,-, a probability of C^* for trading 
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with q Cz Vj via p d Vj, where p might serve as a middle- 
man. Again CI only counts the total number of nodes, 
e.g., triangles, that a node v ^ Vi is connected to, whereas 
Cl* corrects the estimation by taking the total GDP of 
each country into account. We obtain further corrections 
away from linear dependency with the majority of coun- 
tries showing a n.s.i. local cross-clustering coefhcient of 
CI* > 0.9 (Fig.[5j3). This result can be understood by the 
fact that the value of C^* is dominated by nodes with large 
node weights (Eq. Q). Furthermore, in terms of trading 
it is plausible that a country with a large GDP, i.e., high 
Wy, tends to trade with a lot of other countries around 
the globe. Therefore, the probability to form a cluster 
including at least one country of strong GDP is compara- 
tively high, but only those clusters significantly affect C^*. 
Clusters including two partners p,q € Vj with small Wp 
and Wq do not significantly contribute to the value of C^* 
as both weights are multiplied under the sum (Eq. Q). 
Furthermore, as we treat a link not only as a single trade 
relation between countries but rather as an approximation 
of a bundle of links connecting distinct actors in trading 
within the countries, it is clear that the properties of the 
whole bundle are hardly affected by the failure of a sin- 
gle trade connection. Since C^* mirrors this property, we 
again consider C^* an improved representation of the un- 
derlying trade system's properties as compared to Cjj . We 
conclude that the EU-non EU cross-boundary trade net- 
work can be considered as being very robust against the 
failure of links in the network as all nodes tend to cluster 
strongly with their neighbours. As in the first application, 
by introducing a properly chosen typical weight uj in the 
definition of Ci*. further refinements of the measure are 
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Conclusion. — In this work, we have used a gen- 
eral method for transforming standard network measures 
into their weighted counterparts in order to develop node- 
splitting invariant cross-network measures. The latter pro- 
vide an efficient and universal tool for analysing the inter- 
action structure between subnetworks in a network of net- 
works with inhomogeneous node weights rendering those 
measures more appropriate to study real- world problems. 

We have proceeded in three steps: (i) We have car- 
ried out the derivations for two local cross-network mea- 
sures, n.s.i. cross-degree and n.s.i. local cross-clustering 
coefficient. It is important to emphasise that all network 
measures (local and global) can be transformed into their 

(ii) 
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Con- 



node-weighted correspondents as well 
sidering a prototypical spatial model network consisting 
of two interdependent subnetworks embedded on surfaces 
parallel to the Earth's surface, we have validated the be- 
haviour of our two measures of choice. Both n.s.i. cross- 
degree as well as n.s.i. local cross-clustering coefficient 
have been shown to provide an improved representation of 
the model system's underlying network topology as com- 
pared to their unweighted analogues, (iii) In order to illus- 
trate the variety of possible applications of our approach. 



we have carried out an analysis of the complex network 
structure of EU - non EU cross-boundary trade. It was 
shown that node-weighted cross-degree yields additional 
insights into the trade relations as it identifies countries 
showing strong economic dependencies on overseas trade 
such as Great Britain as well as countries driving the econ- 
omy of the EU by exporting to countries with a strong 
economy such as Brazil. In addition, we find that the 
n.s.i. local cross-clustering coefficient indicates a high ro- 
bustness of the trade network to partial failure of links in 
the network, e.g., the breaking of single trade relationships 
between actors in different countries. 

Both applications point out that the n.s.i. cross-network 
measures provide a general tool for analysing networks 
of networks with heterogeneously distributed node impor- 
tance, size or weight that appear in research fields as di- 
verse as climatology, neuroscience or economics. These 
network measures enable us to provide an improved un- 
derstanding of the underlying complex system as long as 
weights are chosen appropriately to refiect the share of the 
domain of interest that is represented by the nodes. 
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